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Abstract 

We present a general framework for the calculation of soft functions for SCETj observables 
through next-to-next-to-leading order (NNLO) in the strong coupling constant. As an example 
of our formalism we show how it can be used to obtain the complete NNLO soft function for 
the A^-jettiness event shape variable. We present numerical results for two examples with phe¬ 
nomenological impact: the one-jettiness soft function for both electron-proton and proton-proton 
collisions. 
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I. INTRODUCTION 


The physics program at the Large Hadron Collider (LHC) and at other experiments 
increasingly relies npon precision calcnlations within the Standard Model (SM) in order 
to search for small deviations indicative of new physics. At the LHC, Rnn I was marked 
by the discovery and initial characterization of the Higgs boson. Run H will focus on the 
detailed investigation of this new state, in which an evermore precise characterization of the 
SM benchmark will be critical. The small errors for most experimental measurements make 
higher-order QCD calculations mandatory in interpreting the data. Such computations may 
be performed at either hxed-order in QCD perturbation theory, or may additionally include 
the resummation of large logarithms in certain regions of phase space, either through analytic 
resummation or via the use of parton-shower simulations. 

One feature of recent progress in precision calculations is the impact of analytic resum¬ 
mation techniques on our understanding of jet properties. They led to the invention of 
new variables such as A^-jettiness [T] and A^-subjettiness [2] that describe jet substructure. 
Additionally, resummation of large logarithmic corrections improved our description of the 
theoretical treatment of Higgs production in exclusive jet bins [SHIS]. The starting point for 
analytic resummation is a factorization theorem describing the observable under considera¬ 
tion, usually in a region of phase space where an expansion of the full QCD result is possible. 
A typical factorized cross section takes on the schematic form 


(7 ~ 




Yl^n 


( 1 ) 


Here, H describes the effect of hard radiation, B encodes the effect of radiation collinear 
to one of the two initial beam directions, S describes the soft radiation, and Jn contains 
the radiation collinear to a final-state jet. Depending on the observable and process under 
consideration, only a subset of these terms may be present. Perturbative corrections to each 
of these functions are minimized by the appropriate renormalization scale choice. Renor¬ 
malization group equations for each separate function allow these scales to be evolved to 
a common one, in the process resumming large logarithms. We have used the language of 
soft-collinear effective theory (SCET) [T6l - [2n] in our description, although similar quantities 
appear in other approaches to resummation. 

Improving the accuracy of resummation requires both knowledge of the anomalous di¬ 
mensions controlling the evolution of the various functions in the factorization theorem, and 
the perturbative expansion of these quantities to higher orders in the strong coupling con¬ 
stant as- Knowledge of the singular structure of QCD gained with resummation formulae 
has also improved our ability to calculate fixed-order QCD quantities to higher precision. A 
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well-known example is the use of gT-subtraction to calculate cross sections through next-to- 
next-to-leading order (NNLO) [ 21 ] • Currently, the hard function H is known to NNLO for 
numerous phenomenologically interesting processes containing hnal-state jets |22H2Z]. The 
NNLO beam functions for observables such as jettiness and beam thrust are available [2S1I2U] . 
as are the NNLO hnal-state jet functions [30113I]. The only missing ingredient of the fac¬ 
torization formula at NNLO for observables such as jettiness and beam thrust is the soft 
function. The knowledge of these soft functions do not only enable the resummation accu¬ 
racy of the corresponding observables to be improved; they are also important components 
of the recently proposed jettiness-subtraction scheme for the calculation of jet cross sections 
through NNLO in QCD [32]. 

We present in this manuscript a general method of computing soft functions through 
NNLO for SCETi observables. It uses sector decomposition [Ml - 136] to extract singularities 
from the integrals which occur in the calculation and reduce them to a form amenable 
to numerical integration. We illustrate our techniques using the Wjettiness event shape 
variable Tn in as an example. We validate our approach using against known results in the 
literature. The iV-jettiness soft function contains the logarithms ln(77v); and a contribution 
of the form 6{Tn)- The logarithmic corrections at NNLO can be obtained by expanding 
the resummed expression for the soft function to We demonstrate that we reproduce 

these known results with our technique. Our computation of the 5{'Tn) correction is new. We 
present numerical results for two selected examples of recent phenomenological relevant: the 
one-jettiness soft function in electron-proton collisions, and the one-jettiness soft function in 
proton-proton collisions. 

Our manuscript is organized as follows. We review the dehnition of the iV-jettiness event 
shape variable in Section |TT| Our calculational framework is presented in Section where 
we show how to reduce the NNLO soft function to a form suitable for numerical evaluation. 


In Section IV we present numerical results for two examples: one-jettiness in electron-proton 
collisions, and one-jettiness in proton-proton collisions. We conclude in Section [V] 


II. DESCRIPTION OF JETTINESS 


We begin with a brief review of the V-jettiness event-shape variable Tn of Ref. m- Tn is 
dehned by 

’ ■ qk 


7iv = ^ mini 


Q^ 


( 2 ) 


Here, the pi are light-like reference vectors for each of the initial beams and hnal-state jets in 
the problem, while the qk denotes the four-momentum of hnal-state radiation radiation. The 
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Qi are dimensionful variables that characterize the hardness of the beam-jets and hnal-state 
jets. For simplicity, we will set Qi = 2Ei, twice the energy of each jet. Writing the jet 
momenta as pi = Eifii, we have 

Tv = ^ mini {rii ■ Qk} ■ (3) 

k 

We will consider the calculation of Tv through NNLO in QCD, which will receive contribu¬ 
tions from single emission and double emission processes. The contributions of single and 
double-real emission processes to iV-jettiness are given explicitly by the following expressions: 

• Single-real emission: Tv = min* {uj ■ qi}; 

• Double-real emission: Tv = min* {uj ■ qi} + min^ {nj ■ ^ 2 }- 

The physical content of the above expression is that jettiness partitions the phase space of 
each emission according to which external direction it is nearest. A pictorial representation 
of this is given in Fig. [T] In each of the regions, Tv is dehned differently in terms of the 
radiation four-momentum. This leads to the insertion of the following measurement functions 
into the phase space for single and double-real emission processes: 

N N 

= >'ta.®) = ^ea, (4) 

i=l hj=^ 

where we have abbreviated 

QI = 5{Tn - Ur ■ qi) 0{nk ■ qi-Ur- qi), 

Qlj = ^(Jn - Hr ■ qi - Us ■ qj) 6{nk -qi-rir- qi) 9{ni ■ qj - Us ■ qj). ( 5 ) 

k^r l^s 

For the QCD processes considered here, ©jj is symmetric under interchange of either its 
upper or lower indices. This allows us to reduce the number of phase-space regions relevant 
for the calculation of double-real emission processes from nine to six. 


III. CALCULATIONAL FRAMEWORK 

With the variable Tv and the structure of the measurement function discussed, we are 
now ready to discuss the calculation of the soft function. The soft function can be expanded 
as a perturbation series in the strong coupling constant, 

S{Tn) = S^^\Tn) + ^S^^KTn) + (6) 
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FIG. 1. Sketch of the division of phase space into regions for the one-jettiness variable, i, j and 
k denote representative hard directions. A two-dimensional projection of the full space has been 
performed for simplicity of presentation. 

We have suppressed the dependence on the renormalization scale /x. The leading-order 
result is just S{Tn), while the calculation of the NLO contribution S^^\Tn) has been 

discussed extensively in Ref. [2Z]. We focus our attention on the computation of S^‘^\Tn)- 
The diagrammatic contributions to the integrand involve emission of gluons from eikonal 
lines, but are most easily obtained from known results for the NNLO soft limits of QCD 
amplitudes. As with all NNLO calculations there are contributions from two-loop virtual 
corrections, one-loop virtual corrections to single-real emission processes (real-virtual), and 
double-real emission diagrams. The two-loop virtual corrections are scaleless in dimensional 
regularization, leaving only the real-virtual and double-real corrections. We are left with the 
following pieces to calculate: 

• the real-virtual corrections to the single-gluon emission process; 

• the qq double-real emission correction; 

• the double-real gluon emission contribution. 


A. The real-virtual correction 

We begin by discussing the real-virtual contribution to the soft function. It receives 
contributions from diagrams of the form shown in Fig. The integrand resulting from these 
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diagrams can be obtained from the one-loop soft-glnon cnrrent in QCD [38]. We write the 
real-virtnal part of the soft fnnction as 





(7) 


where the explicit form of the integrand E^^{qi) can be obtained from Eq. (26) of Ref. [38] . 
Althongh the general strnctnre of this expression is complex, the strnctnre of QCD at NNLO 
gnarantees that the resnlt takes the form of a snm of emissions of qi from a dipole pair 
(i, j), where z, j denote two hard directions in the problem, together with appropriate color 
correlations. We are therefore led to consider the integration of the following bnilding blocks 
from which the real-virtnal corrections for Tn can be constrncted: 


Sijiqi) 




T\1 - e)T^{l + e) 

r2(l-2e)r(l + 2e) 




rzj • Uj 

2ni-qi Uj ■ qi' 


( 8 ) 


Ca = 3 is the nsnal QCD color constant. 7 ^ is the Enler constant, which arises from 
rewriting the bare conpling constant in terms of the renormalized one. Since the bnilding 
block for the soft fnnction is the integral of the Xij{qi) over the real-emission phase space, 
we will consider the anxiliary qnantity 


r(2),0 

'■RV 


{Tn)^ 


Jd—1, 




(9) 


from which we can form the entire integrated real-virtnal correction. 




j 


N ^ ^ ~ 



N 



FIG. 2. Representative diagrams contributing to the real-virtual piece of the soft function. The 
dashed lines represent the eikonal directions. 
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We now discuss the appropriate representation of phase space for the calculation of this 
integral. We introduce a Sudakov decomposition of the radiation momentum qi in terms of 
two light-like directions rim and n^: 

+ qi^^ + q^i^ (10) 

with = Un ■ qi, qi = rim ■ qi, and rim ■ qi± = ■ qi± = 0. With this decomposition it is 

straightforward to write the following expression for the phase space: 


r ^qi _ VLd- 2 , 1 ( ^m • 2-26 

J (27r)'^-i 4(27r)‘^-i V 2 / ^ 


d^dsd(j)i^^ ^■'■‘^sin (11) 


where we have set q~^ = Tn^, q~ = Tn^/s, and have used the angle to parameterize the 
orientation of qi± in the azimuthal plane. 

Upon plugging in the measurement function of Eq. (|^ and the integrand Iij{qi) into 
the phase space, we arrive at three terms, depending on which 0^ occurs. These form two 
distinct sets: two integrals in which r is one of the two directions i,j appearing in the 
integrand Xij, and one in which it does not. We consider the two integrals corresponding to 
r = i and r = k with k distinct from i,j. The r = i case can be written as 



d‘^ ^qi 
(27r)'^“^ 


01. 


( 12 ) 


where we have introduced the subscript i to denote this contribution. The case r = j can be 
obtained by simply permuting the indices i and j in this result. It is convenient to choose 
the light-cone directions rim = rin = rij for this integral. Doing so, it is straightforward 
to derive the following hnal expression for the integral: 

ktlbrA.) = 2-‘-“ 

X JJ 6* - s]. 


dsdx 2 S ^+^''sin 


(</>i) 


(13) 


We have set 0i = 27ra;2, and have introduced the angles (j)ik that denote the separation 
between qi and the hard directions k in the transverse plane. We have also introduced the 
quantities 


^ 27r^ 14 . n vr^ 4 

Brv = 1-+ —e , 

3 3 15 

, . Ui-rik rij-rik 

Aij^k[x,4>) = -——I -2cos(0)- 


rii ■ rij 


rii ■ rij 


rii ■ rij 


(14) 
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Poles in e occur in three places in from the explicit overall factor of 1/e^; from the 

term in the limit Tn —t 0; from the limit s —)■ 0 of the term in the integrand. 

The presence of the theta function makes this expression difficult to integrate analytically 
and extract the pole in s. However, the poles can be easily extracted using plus-distribution 
expansions in s and Tn- 


X 


-l+e 


1 __ 

= t'5W + T-t 


n=0 


In'^x 


X 


(15) 


where x denotes either s or Tn- One this is done the resulting coefficients of the Laurent 
expansion in e can be easily integrated numerically. The case r = j is easily obtained by 
permuting the indices i and j in the quantity Aij^k- We note that the s —)■ 0 limit is associated 
with an ultraviolet singularity, since q~ ^ oo. In SCETi such singularities can always be 
regulated in dimensional regularization, mapped to the unit hypercube and extracted with 
a variable change of the type used here. 

The second case with r = k proceeds similarly, except that the Sudakov decomposition of 
the radiation momentum instead uses rim = Uk, Un = Ui. We denote this contribution with 
the subscript k. Proceeding as before, we derive the hnal result 


r(2),P 

^RV,k 


(Tn) 


- /N 


{ui ■ {ui ■ Uk) 


— l+e r) —1—4e 


ds ( 1 x 2 3^" sin 


(0l) 


X [Akij{s,(j)i)] ^ - s]. 

l^i,k 


( 16 ) 


This form is again suitable for numerical implementation, as the quantity Akij remains 
hnite throughout the allowed phase space. Using and I^^yk-, the entire real-virtual 

contribution to the NNLO soft function for iV-jettiness can be derived. 


B. The qq double-real correction 


We next consider the correction arising from the emission of a qq pair from the hard 
eikonal lines. Several representative diagrams are shown in Fig. The integrand can again 
be obtained from the QCD result for the emission of a soft qq pair, as presented in Ref. [39] . 
It is expressed in terms of the function 

Pi ■ qiPj ■ q2 + Pj ■ qip q2 - Pi ■ qjqi ■ q2 


Tj{qi,q2) = 


(17) 


(gi ■ q2y[Pi ■ {qi + q2)] [pj ■ {qi + ^2)] 

However, color conservation of QCD amplitudes restricts the ways in which enters the 
soft function. The factorization of the QCD amplitude in the double-soft limit indicates that 


\M{. .. , qi, q2)\‘^ ~ {M\ 


Tjj Tj ■ 


\M), 


(18) 






where we have used color-space notation m in writing the color-correlated product of am¬ 
plitudes on the right-hand side. Color conservation allows us to write 


ETilV<) = 0- 


(19) 


We use this relation to remove all color structures of the form Tj ■ Tj in Eq. (18), by dotting 


Eq. (19) with Tj and solving for the Tj ■ Tj term. Doing so, we find that the coefficient of 


each remaining color structure Tj ■ Tj, with i ^ j, contains the combination 


cTjjr Tjj 


-JJ 


21 , 


IJ- 


( 20 ) 


This is the basic building block of the qq contribution whose integration over phase space 
we will study. 



FIG. 3. Representative diagrams arising from emission of a qq pair. The dashed lines represent the 
eikonal directions. 


The integral we must consider is 





d'^ ^qi d'^ ^q2 
{27iY~^ {27iY~^ 


Jij{(li,(l2)M{qi,q2). 


( 21 ) 


The overall numerical factor comes from expressing the bare coupling constant in terms of 
the renoramlized one. We have kept explicit the overall coefficient of (as/(27r))^ in order 
to make clear the normalization of our result. It is convenient to divide this integral into 
several different structures, according to which phase-space parameterization is most suitable 
for performing the extraction of singularities. We hrst divide it into regions according to 
whether the radiated quarks qi and q 2 are closest to one of the emitting eikonal lines i,j, or 
are closer to non-emitting lines which we label as k, 1. This leads us to five distinct regions 
to investigate: 
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1 . both qi and q 2 closest to the same emitting direction i, which we denote as 

2 . qi and q 2 closest to different emitting directions, which we denote as Iqff/] 

3. qi closest to i and q 2 nearest to k, which we call 

4. both qi and q 2 closest to A;, which we call 

5. qi and q 2 closest to different non-emitting directions, which we denote as ■ 

We furthermore find it convenient to divide the integrand Jij into two structures according 
to whether the denominator is quadratic or linear in the invariant qi ■ q 2 - We label these as 
/ and //, respectively, so that . Written explicitly, the integrands for these 

two structures are 


rri ^ g [Pi ■ QiPj ■ Q2 + Pj ■ qiP-q2? 

(gi ■ q2y[pi ■ {qi + q2W[Pj ■ (gi + 52 )]^’ 

jii ^ 2_ Pt • Pj _ 

’’ (91 ■ 92)|Pi ■ (91 + 92)||Pi ■ (91 + 92)1' 


( 22 ) 


This leaves us with a total of ten integrals to compute. The entire qq contribution to the soft 
function can be obtained by appropriately permuting the indices of these structures. Since 
the computation proceeds similarly for all ten terms, we will focus on the representative 
example which exhibits all of the complexities that must be addressed in the general 

case. 


We begin by performing a Sudakov decomposition of both qi and q 2 as in Eq. (10), 
choosing the light-cone directions m = k and n = i. The explicit representations of the 
transverse vectors qi± and q 2 ± are as follows: 


qi± = \qi±\ (cos(0i),sin(0i);O) , 

= l?2±| (cos((;/)2),sin(02)cos(a);sin((/)2)sin(Q;)hJ . (23) 


The last component of each momentum appearing after the semi-colon denotes the e- 
dimensional component of transverse momentum (We recall that the dimensionality of the 
transverse plane is 2 — 2e in dimensional regularization), denotes a unit vector in the 
—2e-dimensional space. A single angle is needed in the general A^-jettiness case to param¬ 
eterize this direction. It is now straightforward to write down the following expression for 
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the integral: 


r(2),*J,r/ „i_8£ 

^qq,kk ^ ^ 


'rr j • nj [rii ■ Uk] ds dt dx^ dx^ dxQ dQ^""^ 


x['^(l-0] lA(l-A)] "sin "(0i)[-ea;6 "](1 - xg) "[st] +"|s - t 


-l- 2 e 


X 


{(v^ - ^/iy + 4Aa/^} 


2 e 


^ tAkij{s, 0i) + (1 - ^)sAkij{t, (j)2k) ^t+{l- ^)s 

X n ^ [-4fci,Ks, 0n) - s] 6' [Aki,i{t, (j)2i) - t ], 


l^i,k 


where we have introduced the abbreviation 


( 24 ) 


„ TT^ 9 8 ^ q Tt"^ 4 

Brr = 1 — —6 — A —e • 

3 3^ 90 


(25) 


We have made the following variable changes to arrive at this expression: 


+ n- c - + n- n c\ - “ 0 

Qi = iNi, Qi = -, qt = Tn{ 1 - 0: q2 = -7-: 

s t 

01 = 27rx4, A = sin^(7rx5/2), cos(q;) = 1 — 2 xq . 


(26) 


We have in addition followed the sector decomposition approach to NNLO calculations 
HD 112] and have made a non-linear change of variables to map cos( 02 ) to the unit hypercube. 


The quantity denotes the angular parameterization of the direction in Eq. (23), 


normalized so that it integrates to unity. We note that in the zero-jettiness and one-jettiness 
cases, we may immediately integrate over to obtain unity. 

This integral is not yet suitable for numerical implementation, as there are singularities 
associated with the joint limit s, f —?• 0 that cannot yet be extracted with a plus-distribution 
expansion. We order these two limits by inserting the following partition into phase space, 
as is done in sector decomposition of real radiation in QCD [HI HU 112] : 


1 = 6{s — t) + 9{t — s). 


(27) 


We then remap the limits of integration to the unit hypercube. Doing so renders all singu¬ 
larities amenable to a plus distribution example. Focusing on the s > t sector for illustrative 
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purposes and setting ^ = xi, s = X 2 , and t = X 2 (l — X 3 ), we find the following final result: 


^ qq,kk,s>t ^ ^ 


RR 


'N ni-nj[ni-nk\ 


dxi dx2 dxz dxi dx^ dx^ 


\l+e-l-2€ 


x[xi(l-a:i)] ^"[A(l-A)] "sin ^"(0i)[-exg^ "](1-X6) "xf (1 - X 3 )^+"x 3 


X 


{(1 - \/l - xsf + 4Aa/1 - X 3 } 


2 e 


Xi(l - X 3 ) tAkij{x2, (j)i) + (1 - xi)Akij{x 2 {l - X 3 ), (j)2k) 1 - a;iX 3 

X JJ 6* [Aki,i{x2, 0 n) - X2] 61 [Afci,,(x2(l - X3), 02 i) - 3:2(1 - 2:3)]. 


l^i,k 


(28) 


The poles in this expression come only from expanding the overall TA 


-l- 4 e 

N 


and the factor 


„-l- 2 e 


Xq 


using the plus distribution expansion of Eq. (15). Once this expansion in performed, 
the integral can be simply evaluating numerically. The computations of the other sector and 
the other integrals required for the qq contribution proceed similarly. Only the ordering of 
the s and t limits is needed to make all singularities manifest. 


C. The gg double-real correction 


Finally, we consider the correction arising from the emission of two gluons from the hard 
eikonal lines. Several representative diagrams are shown in Fig. The integrand can be 
obtained from Ref. [31] • Two distinct types of contributions to the soft function can be 
identified. The first type comes from the squares of single-gluon currents, and can be easily 
obtained following the techniques of Ref. [3Z]- We do not discuss it further here. The second 
type are genuinely non-abelian contributions proportional to Ca that are not simply the 


square of one-loop terms. These can be written in the form of Eq. (|18|), except with the 

The 


replacement of eikonal functions lij —?■ iS^-, with Sij given in Eq. (110) of Ref. 
same color conservation argument as for the qq case indicates that only the combination 


l~ij — Sii + Si-i — 2 S, 


jj 


(29) 


contributes to the result. Simple algebraic manipulation leads us to the result 

Tij = (1 - e)j/j + 2j// + 


q(s.o.) _ Pi ' Pj 
^ij ~ 


+ 


Pi ■ qiPj ■ q2 + Pj ■ qiPi ■ q2 A 
[P*-(?i + ?2)][p,-(gi + g2)] J ’ 
1 ^ {Pi- Pj? 


qi ■ q2 \Pi ■ qiPj ■ q2 Pj ■ qiPi ■ q2 


Pi ■ qiPi ■ q2Pj ■ qiPj ■ q2 


(30) 
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The first two terms with and are the same structures as found in the qq case. The 
structure proportional to the function arises from the strongly-ordered limit of QCD, 

in which there is a hierarchy between the energies of the radiated gluons. It is straightforward 


to follow the same parameterizations and steps presented in Sec. Ill B to render this structure 


suitable for numerical evaluation, and we do not repeat the details here. 


"s / 

^ v\/\AA/ 

/ 



' / 


FIG. 4. Representative diagrams arising from emission of a gg pair. The dashed lines represent the 
eikonal directions. 


IV. NUMERICAL RESULTS FOR ONE-JETTINESS 

We present in this section numerical results for the NNLO contributions to iV-jettiness, 
including validation against known results in the literature. We focus on two example cases in 
order to illustrate our results. We begin with one-jettiness in electron-proton collisions, which 
has received recent interest in the contexts of probing nuclear dynamics in electron-nucleus 
collisions |l3l 01], and of improving jet phenomenology in deep inelastic scattering jH] 06 ] . 
We also consider one-jettiness in proton-proton collisions, for which the NNLO soft function 
is a necessary component of a recently introduced subtraction scheme for NNLO fixed-order 
calculations jaa |33]. The restriction to one-jettiness only simplihes the color structure of 
the integrands. The basic building blocks are those presented in the previous section. 


A. One-jettiness in ep collisions 

We begin with a presentation of the one-jettiness soft function in ep collisions as a valida¬ 
tion of our calculation. In this case the soft function can be obtained analytically from the 
known result for the thrust distribution in collisions m- We present the analytic result 
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in the Appendix. It is most convenient to express the result in terms of the soft contribution 
to the cumulative cross section for one-jettiness: 



with L = log(7^™*//i). Since all components of the soft function contain the overall depen¬ 
dence the integration over 7i to obtain this cumulant is trivial to perform. We use 

the numerical approach described in the previous section and compare it against the ana¬ 
lytic result for the Cq coefficient. We compare in Figs. |^[^and[^the separate contributions 
from the coupling-constant renormalization, the double-real emission contribution and the 
real-virtual correction. In the former two cases, we further separate the NpT^ and C^Ca 
color structures. Each contribution is plotted as a function of S 12 = ni ■ 77 - 2 . In all cases, the 
numerical prediction for Cq agrees perfectly with the analytic calculation. 


NfTf^ term from Renormalization term from Renormalization 




FIG. 5. Comparison between the analytic and numerical calculations of the Cq coefficients from 
renormalization, for both the Np Tp (left panel) and Cp Ca (right panel) terms, as a function of 
S 12 . The blue solid lines represent the analytic calculation and the red dots are from our numerical 
approach. 


B. One-jettiness in pp collisions 


We next consider the one-jettiness soft function in proton-proton collisions. We calculate 


the 0{al) cumulative cross section defined in Eq. (31). The one-jettiness soft function 
depends on three hard directions, which we label as ui, 77 , 2 , and 773 . We align 77i and n 2 
with the incoming beam axes in the ±z directions, and let 773 lie along the outgoing jet 


direction. The Ci in Eq. (31) can then be written as functions of the kinematic invariant 
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NfT/^ term from Double-Real 



CfCfy term from Double-Real 



FIG. 6. Comparison between the analytic and numerical calculations of the Co coefficients from 
the double-real contribution, for both the NpTji (left panel) and CfCa (right panel) terms, as a 
function of Si 2 . The blue solid lines represent the analytic predictions and the red dots are from 
our numerical approach. 


C/rC^ term from Real-Virtual 



FIG. 7. Comparison between the analytic and numerical calculations of the Cq coefficients from 
the real-virtual contribution. The blue solid lines represent the analytic predictions and the red 
dots are from our numerical approach. 


Si 3 = ni-rio. We note that G 4 , C 3 , C 2 and Gi can be determined by expanding the resummed 
soft function to order a^. We use these terms to verify our direct hxed-order calculation. We 
show the results for the Abelian and the non-Abelian contributions separately. In both cases 
we assume that the color factors associated with partons in the ni, n 2 and 77,3 directions are 
Cp, Ca and Cp, respectively. 

In Figs. and we compare the coefficients Cn of L" for n = 4, 3, 2,1 from the direct 
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NNLO computation and the expanded resummed soft function to order aj for the Abelian 
and non-Abelian contributions. We hnd complete agreement between the hxed-order calcu¬ 
lations and the resummation predictions. In £g. 1^ we show our calculations of Co for both 
the Abelian and non-Abelian terms. We note that the Abelian result can be predicted from 


exponentation of soft emissions. This prediction is shown in Fig IT, and agrees perfectly 
with our numerical derivation. The non-Abelian contribution to the Cq term is new. 



■S'13 



•J|3 



FIG. 8. Comparison between the analytic and numerical calculation of the coefficients of 
(n = 4,3)2,1) for the Abelian contribution to the soft function. The red solid lines represent our 
direct NNLO calculation, and the blue dots are predicted by expanding the resummed results to 

0(a2). 


V. CONCLUSIONS 

In this manuscript we have presented a calculational framework for obtaining NNLO 
results for the soft functions which are ubiquitous in effective held theory descriptions of 
scattering processes. Using the approach of sector decomposition, which has long been a 
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FIG. 9. Comparison between the analytic and numerical calculation of the coefficients of 
(n = 4, 3, 2,1) for the non-Abelian contribution to the soft function. The red solid lines represent 
our direct NNLO calculation, and the blue dots are predicted by expanding the resummed results 
to 0{al). 


crucial part of the fixed-order QCD calculational toolbox, we have extracted all singular¬ 
ities from the soft-function integrands and reduced the computation to a set of numerical 
integrals which are simple to evaluate. We have illustrated our technique using the N- 
jettiness event shape variable [1]. Several recent ideas suggested in the literature require 
the NNLO soft function for A^-jettiness. These include the improved theoretical description 
of electron-nucleus collisions |13] and deep inelastic scattering aa. and a recently proposed 
NNLO subtraction scheme for LHC processes containing hnal-state jets |32]- We have shown 
numerical results for the one-jettiness soft function in both ep and pp collisions that address 
these needs. 

Although we have focused on iV-jettiness as an example, the techniques introduced hold 
more generally for SCETj observables which feature constraints on the light-cone momenta 
of the hnal-state radiation. In such cases all singularities that appear in the integrand can 
be regulated with dimensional regularization. For SCETn observables which instead have a 
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FIG. 10. contributions to the 0(ag) soft function, for both Abelian (left panel) and non-Abelian 
(right panel) cases. The red solid lines represent our direct NNLO calculation. In the Abelian case, 
the blue dots are obtained by expanding the exponentiated NLO soft function. 

constraint on the transverse momentnm of the measured radiation, additional singularities 
not controlled by dimensional regularization appear. We believe that our formalism can be 
extended to also provide NNLO calculations in such cases when an appropriate regulator is 
chosen [181 - 150] . We look forward to pursuing these and other extensions of our work. 
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APPENDIX 

We present here the analytic result for the one-jettiness soft function in electron-proton 
collisions. We organize this result according to color structure, and to whether it arises from 
a double-real or a real-virtual contribution: 

sW( 7 :) = + (32) 
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The individual components are as follows: 
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We have set Lij = log(nj • nj/2). 
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